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ABSTRACT 

We report on an efficient method for locating the apparent horizon in numerically 
constructed dynamical 3D black hole spacetimes. Instead of solving the zero ex- 
pansion partial differential equation, our method uses a minimization procedure. 
Converting the PDE problem to minimization eliminates the difficulty of imposing 
suitable boundary conditions for the PDE. We demonstrate the effectiveness of this 
method in both 2D and 3D cases. The method is also highly parallelizable for 
implementation in massively parallel computers. 

o 

^ 1 

The apparent horizon (AH) is a crucial characteristic of a black hole spacetime. 
CN For numerically constructed black hole spacetimes, while there are efficient methods 

to find the AH in a wide variety of 2 dimensional situations/ there is as yet no 
^ satisfactory method for the 3D cases. 

Q The by now standard method of finding AH, at least for lower dimensional cases, 

C^l is to solve the zero expansion equation: 

^ = 3»|_J^^ + J^;.^.3V=O , (1) 
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where is the outward spatial normal vector of the AH on the constant time slice at 



O 

. . 

I which the AH is to be found, and A is the trace of the extrinsic curvature A,-, of that 

toX) slice. This equation can be easily integrated to determine the AH in spacetimes with 

symmetries, as the symmetries often provide boundary conditions at the edges of the 
computational domain (e.g., the axis and the equator). However, for a general 3D 
spacetime, as the AH is topologically a 2-sphere, one immediately faces the difficulty 
of imposing suitable boundary conditions even for the starting of the integration, let 
alone the devising of an efficient method for solving this nonlinear partial differential 
equation (PDE). 

Our method for handling this problem of finding the AH surface in a general 3D 
spacetime is to convert the PDE problem to a minimization problem, the solution of 
which is well understood. This can be achieved by: 

1. Express a surface F[xi) = 0, which is topologically a 2-sphere living on a 
constant time slice, in terms of the symmetric trace free (STE) tensors,^ or 
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other suitable complete set of basis functions, denoted by P„: 

N 

F{xi) = J2A.,P., = . (2) 

n=0 

As our 3D codes for evolving the black hole spacetimes are written in cartesian 
coordinates, STF tensors are most convenient for our purpose. For relatively 
smooth surfaces, as one expects the AH to be, only the first few lower rank 
tensors (i.e., N being finite and not too large, see below) are needed. 

2. A surface is hence represented by the set of coefficients A„. The expansion at 
each point of the surface is a function of A„. The condition for the trial surface 
to be the apparent horizon becomes 

/ 0'(A„) = (3) 

As 0^(A„) is semi-positive definite, there are efficient optimization algorithms 
to search the A^-space for the surface closest to the AH among all test surfaces 
so parameterized. 

2. Implementation and Results 

We first applied our AH finder based on this method to 2D testbeds. We used 
as our background spacetime data obtained from a code developed by Bernstein et 
al.^ This code evolves a black hole distorted by an axisymmetric distribution of 
gravitational waves (Brill wave). The black holes can be highly distorted by the 
incoming wave. We compare the results obtained with our new AH finder with the 
results from the AH finders constructed using the standard PDE method. For test 
surfaces with 16 coefficients (with the spherical harmonics Yi^ as basis functions), we 
find that both the new method and the old methods coincide within the given accuracy 
of the PDF solvers. The new method takes less than 5 minimization iterations to 
converge to the correct surface in typical situations where we use the previous solution 
as a guess on the present time slice. 

For the 3D testbeds, we use the symetric trace free (STF) tensors as basis func- 
tions. For the Schwarzschild spacetimes constructed in our 3D cartesian codes, we 
find that out method can converge to the correct Schwarzschild horizon location 
within a few iterations and using only a 6 coefficient expansion. The horizon is cor- 
rectly located to within a small fraction of a grid zone. Fig. 1 shows how the 3D AH 
finder converges. We are presently testing the AH finder on evolved Schwarzschild 
black holes and also on two black hole data sets in 3D cartesian coordinates. Results 
of this work will be reported elsewhere. 

3. Conclusion 

We have developed a promising method of finding the AH in a numerically con- 
structed spacetime by a minimization procedure. This method should be an efficient 
and accurate tool for use in 3D numerical relativity. 




Figure 1: (A) Initial trial surface. (B) After 4 minimization iterations, we converge to the spherical 
Apparent Horizon for a Schwarzschild black hole in a 3D cartesian grid. 
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